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We present the results of a Bayesian analysis of solar and KamLAND 
neutrino data in the framework of three-neutrino mixing. We adopt 
Ph' two approaches for the prior probability distribution of the oscillation 

parameters Am^, sin 2 $12, sin 2 $13: 1) a traditional flat uninformative 
prior; 2) an informative prior which describes the limits on sin 2 $13 
obtained in atmospheric and long-baseline accelerator and reactor neu- 
^ \ trino experiments. In both approaches, we present the allowed regions 

in the sin 2 di^-Am^ and sin 2 f^-sin 2 1? 13 planes, as well as the marginal 
posterior probability distribution of sin 2 $13. We confirm the 1.2<r hint 
■^j- ■ of $13 > found in hep-ph/0806 . 2649 from the analysis of solar and 

KamLAND neutrino data. We found that the statistical significance of 
the hint is reduced to about 0.8<r by the constraints on sin 2 -$ 13 coming 
00 ; from atmospheric and long-baseline accelerator and reactor neutrino 

data, in agreement with hep-ph/0808 . 2016. 

> : 

X ■ 1 Introduction 



Neutrino oscillation experiments have shown that neutrinos are massive and mixed par- 
ticles (see Refs. [1-3]). Solar and KamLAND neutrino experiments observed { v e — > v^ r 
transitions due to neutrino oscillations generated by a squared-mass difference 

Am 2 SOL ~ 8 x 10~ 5 eV 2 . (1.1) 

Atmospheric and long-baseline accelerator neutrino experiments measured'^ — ^v T tran- 
sitions due to neutrino oscillations generated by a squared-mass difference 

Am 2 ATM ~ 2.5 x 10" 3 eV 2 . (1.2) 
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Figure 1: The two three-neutrino schemes allowed by the hierarchy Amg OL <C Am\ TM . 

Hence, there is a hierarchy of squared-mass differences: 

Am ATM - 30 Am^ OL . (1.3) 

This hierarchy is easily accommodated in the framework of three-neutrino mixing, in 
which there are two independent squared-mass differences. We label the neutrino masses 
in order to have 

A^sol = A«4 , (1-4) 

A m ATM - l Am 3ll - \ Am 32\ ■ (1-5) 

The two possible schemes are illustrated in Fig. [TJ They differ by the sign of Am^ ~ 
Am 2 32 . 

For the 3x3 unitary mixing matrix of neutrinos we adopt the standard parameteri- 
zation in Eq. (I A. 21) of Appendix |X] [4,5]. The negative results of the Chooz [6] and Palo 
Verde [7] long-baseline neutrino oscillation experiments, together with the evidence of 
neutrino oscillations in atmospheric and long-baseline accelerator neutrino experiments, 
imply that the mixing angle #13 is small [8] (see Ref. [9] for updated bounds). On the 
other hand, the values of the other two mixing angles are known to be large from the 
results of solar and KamLAND experiments ($12) arid the results of atmospheric and 
long-baseline accelerator neutrino experiments ($23) • 

In Ref. [10] we presented the results of a Bayesian analysis of the solar and Kam- 
LAND neutrino data in the framework of two-neutrino mixing, which is obtained from 
three-neutrino mixing in the approximation of negligible $13. In this paper, we extend 
our Bayesian analysis to the framework of three-neutrino mixing, aiming at the determi- 
nation of the constraints on the value of the small mixing angle $13 implied by solar and 
KamLAND neutrino data. 

The plan of the paper is as follows. In Section [2] we present the constraints on the 
value of $13 in a standard x 2 analysis, to be compared with the Bayesian results with a 
uninformative prior presented in Section [3j In Section H] we present the results obtained 
with an informative prior which represents information on $ 13 obtained in atmospheric 
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Figure 2: The 90%, 95%, and 99.73% C.L. regions in the sin 2 plane obtained 

in the least-squares analysis. The solid and dotted lines enclose, respectively, the regions 
obtained with solar data and KamLAND data. The shadowed areas are obtained from the 
combined analysis of solar and KamLAND data. The figure on the right is an enlargement 
of the interesting area of the figure on the left. The dot, cross, and asterisk indicate, 
respectively, the best-fit points of the solar, KamLAND, and combined analyses. 

and long-baseline accelerator and reactor neutrino experiments, independently from solar 
and KamLAND neutrino data. The conclusions are given in Section El 

2 x 2 Analysis 

The traditional way to extract information on the neutrino mixing parameters from solar 
neutrino data is based on the standard least-squares method, also called "x 2 analysis". 
The least-squares function % 2 is given by x 2 — ~ 2 In £+ const ant, where £ is the likelihood 
function. In this section we present our results in a traditional least-squares analysis of 
solar and KamLAND neutrino data. The least-squares function described in this Section 
will be used in the Bayesian analysis presented in Section [3] for the calculation of the 
sampling probability distribution, which is proportional to the likelihood function. 

We consider the data of the following solar neutrino experiments: Homestake [11], 
GALLEX/GNO [12], SAGE [13], Super-Kamiokande [14] [15], and SNO [16]. The least- 
squares function of solar neutrino data is given by 

xl = E (^ xp - *f ) MT% (£f p - Rf) ■ (2.1) 

Here R* xp are the solar data points, whose number is N$ = 80, accounted as follows: 

• the rate of the Homestake 37 CI experiment [11]; 

• the combined rate of the 71 Ga experiments GALLEX/GNO [12] and SAGE [13], 
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Figure 3: The 90%, 95%, and 99.73% C.L. regions in the sin 2 $i 2 -sin 2 plane obtained 
in the least-squares analysis. The solid and dotted lines enclose, respectively, the regions 
obtained with solar data and KamLAND data. The shadowed areas are obtained from 
the combined analysis of solar and KamLAND data. The dot, cross, and asterisk indicate, 
respectively, the best-fit points of the solar, KamLAND, and combined analyses. 

• the day and night energy spectra of the Super-Kamiokande experiment [14] (21 + 21 
bins) and Super-Kamiokande experiment [15] ; 

• the day and night energy spectra of charged-current events in the SNO experiment 
[16] (17 + 17 bins); 

• the neutral-current event rate in the salt phase of the SNO experiment [16]. 

• the NCD neutral-current event rate in the SNO experiment [17]. 

The corresponding theoretical rates Rf 1 depend on the neutrino oscillation parameters. 
The covariance error matrix Vs takes into account the correlations of theoretical uncer- 
tainties, according to the discussions in Refs. [18-21]. In our analysis, the initial flux of 
8 B solar neutrinos is considered as a free parameter to be determined by the fit, mainly 
through the SNO neutral-current data. For the other solar neutrino fluxes, we assume 
the BP04 Standard Solar Model [22]. The transition probability in the Sun is calculated 
using the standard method [23] based on the hierarchy of squared-mass differences in 
Eq. (11.31) . which implies that the oscillations generated by the large mass-squared differ- 
ence Am\ TM are averaged out (see Refs. [1-3]). For the calculation of the regeneration 
of solar i/ e 's in the Earth, we use Eq. flA. 12[) . derived in Appendix IA1 

Neutrino oscillations due to the same mixing parameters which generate the oscilla- 
tions of solar neutrinos have been observed in the KamLAND very-long-baseline reactor 
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Figure 4: A% 2 as a function of sin 2 $i 3 . The straight horizontal lines show the levels 
corresponding to 90%, 95%, and 99.73% C.L.. 

neutrino oscillation experiment [24]. The KamLAND least-squares function is3 

N K 

Xl = £ (AT P - N?) (V£% (iVf p - Af ) , (2.2) 

where iVx = 17 is the number of energy bins, A^ exp is the number of events measured in 
the zth bin and N* h is the corresponding theoretical value, which depends on the neutrino 
oscillation parameters. The covariance error matrix Vk takes into account the statistical 
uncertainties and the correlated and uncorrelated systematic uncertainties, all added in 
quadrature. 

The global least-squares function is 

Xt = Xs+Xk- (2-3) 

We minimized Xt with respect to the three mixing parameters Am^, sin 2 $12, and 
sin 2 i?i 3 . We found the best-fit point 

Am 2 x = 7.58 x 1(T 5 eV 2 , sin 2 12 = 0.31 , sin 2 13 = 0.021 . (2.4) 

The 90%, 95%, and 99.73% C.L. regions in the sin 2 ^ 13 -Am 2 1 and sin 2 i? 12 -sin 2 1? 13 planes 
are shown, respectively, in Figs. [5] and El These regions correspond to 2 degrees of 

1 In Ref. [10] and in the first version of this paper (arXiv : 0810 . 5443vl) we adopted a different 
least-squares function, which is appropriate for a Poisson distribution (see Refs. [5,25]). We think that 
the Gaussian least-squares function in Eq. (|2.2[) is more appropriate for the analysis of KamLAND data, 
because it allows us to take into account the systematic uncertainty in each energy bin, as discussed in 
Ref. [26]. 
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Figure 5: The 90%, 95%, and 99.73% Bayesian credible regions in the sin 2 
plane obtained with an uninformative constant prior probability distribution. The solid 
and dotted lines enclose, respectively, the credible regions obtained with solar data and 
KamLAND data. The shadowed obtained from the combined analysis of solar 

and KamLAND data. The figure on the right is an enlargement of the interesting area 
of the figure on the left. The dot, cross, and asterisk indicate, respectively, the best-fit 
points of the solar, KamLAND, and combined analyses. 

freedom. The third parameter ($ 12 in Fig. [5] and Am 2 2 in Fig. [3]) is marginalized by 
minimizing Xt- 

From Figs. [2] and El one can see that KamLAND data constrain $13 more than solar 
data. 

Figure H] shows the difference A% 2 of x 2 from its minimum value as a function of 
sin 2 $13. The resulting 90%, 95%, and 99.73% C.L. upper bounds for sin 2 $ 13 , determined 
by the intersection of the A% 2 curve with the straight horizontal lines in Fig. HI are, 
respectively, 



It is interesting to note that the best-fit point for sin 2 $13 in Eq. (12.4D is slightly 
larger than zero, in agreement with the value obtained in Ref. [27] (see also Refs. [9,28]), 
sin 2 $ 13 = 0.021 ± 0.017. Since we have sin 2 tf 13 = 0.021 ± 0.018, our hint of i7 13 > is at 
the 1.2cr level (the precise value of A% 2 for -$ 13 = is 1.33, corresponding to 1.15a"). 

3 Bayesian Analysis 

In the Bayesian approach, the analysis of the data allows us to calculate the posterior 
probability distribution of the mixing parameters, assuming a prior probability distribu- 
tion which quantifies the prior knowledge. Denoting with M. = {Am 2 !, sin 2 1? 12 , sin 2 -$ 13 } 
the set of mixing parameters to be determined by the analysis, the normalized posterior 



sin 2 i7i3 < 0.051 (90%) , 0.057 (95%) , 0.076 (99.73%) . 



(2.5) 
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Figure 6: The 90%, 95%, and 99.73% Bayesian credible regions in the sin 2 t?i 2 -sin 2 i?i 3 
plane obtained with an uninformative constant prior probability distribution. The solid 
and dotted lines enclose, respectively, the regions obtained with solar data and Kam- 
LAND data. The shadowed areas are obtained from the combined analysis of solar and 
KamLAND data. The dot, cross, and asterisk indicate, respectively, the best-fit points 
of the solar, KamLAND, and combined analyses. 



probability distribution of the mixing parameters is given by 

P(M\V,1)= IdMp{vlM ^ )p{Mliy (3.1) 

where p(T>\Ai,X) is the sampling probability distribution, p(Ai\I) is the prior probability 
distribution, and &M. = dAm^ d sin 2 $12 dsin 2 $13. The symbols T> and X represent, 
respectively, the experimental data and all the prior general knowledge and assumptions 
on solar and neutrino physics. 

The sampling probability distribution is given by 

p(V\M,I) oc (\V s \\V K \y 1/2 e~^ 2 , (3.2) 

with the least-squares function Xt given in Eq. (I2.3p . The normalization factor is irrele- 
vant, since it cancels in Eq. ( 13. ip . We retained only the coefficient (|Vjg| |^k|) -1 ' 2 , which 
depends on the neutrino mixing parameters in Ai (see Ref. [21]). 

In Ref. [10] we have shown that the choices of constant uninformative priors in the 
sin 2 ^-Am 2 ,! or log sin 2 $i 2 -log Am 2 x planes are practically equivalent, because of the 
excellent quality of the data. Hence, in the three-neutrino mixing analysis we assume 
a constant prior in the three-dimensional space of the parameters Am^, sin 2 $12, and 
sin 2 $ 13 . 
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Figure 7: Comparison of the allowed regions in the sin 2 i^-Am 2 ^ plane obtained with the 
X 2 analysis and the Bayesian approach. The light- shadowed and light+dark-shadowed 
areas cover, respectively, the Bayesian credible regions with £ probability and the x 2 
region with £ C.L.. In the three figures, from left to right, £ = 90%, 95%, 99.73%. 



Figures [5] and [6] show the resulting credible regions with 90%, 95%, and 99.73% prob- 
ability in the sin 2 ^-Am 2 ^ and sin 2 $i 2 -sin 2 $13 planes, respectively. A credible region is 
the smallest region with the given integral posterior probability. In practice, a credible re- 
gion is calculated as the two-dimensional region surrounded by an isoprobability contour 
which contains the point of highest posterior probability. In each plane of parameters, the 
probability distribution is calculated by integrating p(Ai\T>,X) over the third parameter 
(sin 2 $12 in the plane sin 2 and Am 2 2 in the plane sin 2 ^12-sin 2 $13). 

The credible regions in Figs. El and El are similar but slightly lager than the x 2 -allowed 
regions in Figs. [5] and [31 The comparison of the two types of region is shown in Fig. Hand 
El It is fair to conclude that the Bayesian analysis with an uninformative prior confirms 
the results obtained with the traditional x 2 analysis. 

Figure [9] shows the marginal posterior probability distribution of sin 2 $13, which is 
given by 

p(sin 2 tf 13 |£>,J) = J dAm 2 21 J dsm 2 $ 12 p(M\V,l) . (3.3) 

The resulting credible upper bounds for sin 2 ^ 13 with 90%, 95%, and 99.73% probability 
are, respectively, 

sin 2 $i3 < 0.048 (90%) , 0.054 (95%) , 0.075 (99.73%) . (3.4) 

These limits are similar to those in Eq. (12.51) . in agreement with the above conclusion 
that an uninformative-prior Bayesian analysis confirms the results obtained with a x 2 
analysis. 

We investigated also the hint of $13 > in the Bayesian approach. Since the proba- 
bility of the smallest posterior credible region which includes $13 = is 0.86, considering 
the rescaled probability corresponding to a two-tailed posterior Gaussian distribution we 
obtain a hint of $13 > at the 1.2a level, as in the x 2 analysis (see the discussion at the 
end of Section |2J) . 

As a caveat on the comparison of frequentist and Bayesian results, let us remind 
that the two theories are based on different definitions of probability. Hence, although 
"numerical results tend to be the same for the two approaches in the asymptotic regime, 
that is, when there are a lot of data, and statistical uncertainties are small compared 
with the distance to the nearest physical boundary" [29], the interpretation is different. 
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Figure 8: Comparison of the allowed regions in the sin 2 ^12-sin 2 $13 plane obtained 
with the x 2 analysis and the Bayesian approach. The light-shadowed and light+dark- 
shadowed areas cover, respectively, the Bayesian credible regions with £ probability and 
the x 2 region with £ C.L.. In the three figures, from left to right, £ = 90%, 95%, 99.73%. 

In the frequentist theory the probability of a class of random events is the relative 
frequency of occurrence of these events when the total number of events tends to infinity. 
In parameter estimation, a confidence interval with a CL. is an element of a hypothetical 
set of confidence intervals which have a frequentist probability a of covering the true value 
of the parameter (see Refs. [29,30]). Notice that in frequentist statistics it is not allowed 
to make any statement about the true value of the parameter, which is a fixed unknown 
number, not a random variable, albeit in practice frequentist statistics is very often 
applied to quantities which are not random variables, as systematic errors. The correct 
frequentist statements in parameter estimation concern intervals in the parameter space 
and the frequency of their coverage of the unknown true value in the asymptotic limit. 
This is the meaning of the allowed regions in Figs. [2] and [3j The 1.2<r hint of $13 > 
discussed at the end of Section [2] means that the best-fit value of sin 2 ^ is 1.2a away 
from sin 2 ?9 1 3 = 0, i.e. the confidence intervals obtained in the x 2 analysis with less than 
about 76% CL. do not include sin 2 -# 13 = 0. 

In the Bayesian theory probability represents the degree of belief based on the available 
knowledge. Hence it is possible to estimate a probability for any kind of event, not only 
for random variables as in frequentist statistics. In particular, systematic errors can be 
treated without any inconsistency. Moreover, there is no need to consider hypothetical 
quantities, since the posterior probability distribution is straightforwardly obtained from 
the prior probability distribution and the sampling probability distribution using Bayes' 
theorem, as in Eq. (13.11) . The only difficult task in Bayesian theory probability is the 
estimation of the prior probability distribution on the basis of the available knowledge. 
In parameter estimation, one can calculate the Bayesian probability of the true value of 
the parameter to lie in an interval by integrating the posterior probability distribution. 
The 1.2a hint of $13 > discussed above means that the credible intervals obtained with 
a Gaussian approximation of the posterior probability density having less than about 
76% probability do not include sin 2 $13 = 0. Note that the Gaussian approximation 
of the posterior probability density is defined on the whole real axis of sin 2 1? 13 for the 
comparison with the analogous frequentist result using the traditional terminology In 
fact, the least-squares analysis leads to correct frequentist confidence intervals only in 
the case of a Gaussian likelihood in which the mean values of the data points are linear 
functions of the parameters. In practice this requirement is approximately satisfied in a 
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Figure 9: Marginal posterior probability distribution of sin 2 $13 obtained with an un- 
informative constant prior probability distribution. The straight vertical lines show the 
levels of 90%, 95%, and 99.73% integrated probability. 

region around the minimum of the x 2 if the data are abundant and the minimum of the 
X 2 lies far from any boundary of the parameters. Since in the case under consideration 
we are close to the boundary sin 2 -^ > 0, we can compare the Bayesian result with the 
frequentist least-squares result in which the boundary has not been taken into account 
only by relaxing the boundary restriction. 

4 An Informative Prior 

In the previous Section we analyzed the solar and KamLAND neutrino data assuming 
a constant uninformative prior probability distribution in the three-dimensional space 
of the parameters Am| 1; sin 2 $12, and sin 2 $13. However, as remarked in the introduc- 
tory Section [TJ the value of $13 was known to be small before the analysis of solar and 
KamLAND neutrino data from the negative results of the Chooz [6] and Palo Verde [7] 
long-baseline neutrino oscillation experiments combined with the evidence of neutrino 
oscillations in atmospheric and long-baseline accelerator neutrino experiments. In the 
Bayesian approach it is natural to try to express this prior knowledge through a prior 
probability distribution. The resulting posterior probability distribution of the mixing 
parameters Am 2 1; sin 2 $ 12 , and sin 2 $13 is interpreted as our knowledge about their values 
obtained from solar and KamLAND neutrino data, taking into account the information 
on $13 obtained in atmospheric and long-baseline accelerator and reactor neutrino exper- 
iments. 

Since we do not have the machinery for the fit of the data of atmospheric and long- 
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baseline accelerator and reactor neutrino experiments, we constructed a prior probability 
distribution for $13 using the x 2 reported in Fig. 24 of Ref. [31], where such fit was 
performed. In Fig. 24 of Ref. [31] there are two slightly different curves corresponding to 
the normal and inverted schemes (see Fig. [T]), which depict x 2 (cos5sin$i3) f° r the t wo 
CP-conserving cases cos 5 = ±1, where 5 is the phase in the mixing matrix in Eq. (IA.2I) . 
Since we do not have any information on the value of 6, for each scheme we considered a 
prior probability distribution for $13 marginalized over cos 5 = ±1: 

/ n n-x ( X 2 (cos5sin$ 13 )\ 
p{$i 3 \T) oc 2^ ex P( 2 ]' ^ ' 

cos<5=±l ^ ' 

For sin 2 $ 12 and Am^ we assumed constant uninformative priors as in Section [HI 

The prior distributions (14. ID in the normal and inverted schemes are depicted by 
the dotted curves in Fig. [12], from which one can see that they have a maximum for 
sin 2 ^^ = 0. Hence, they disfavor the hint of $13 > 0. The 90%, 95% and 99.73% prior 
upper bounds for sin 2 $ 13 are 

• 2 / 0.030 , 0.036 , 0.051 (normal scheme) , , , 

sin t> 13 < I Q Q33 ^ Q 03g ^ Q Q57 (inverted scheme) . ^ ' ' 

Notice that such disfavoring of the hint of $13 > obtained from the x 2 i n Fig- 24 of 
Ref. [31] in the Bayesian approach is in contrast with a faint hint of $i 3 > which can be 
obtained in the frequentist approach by considering the minimum of x 2 a t cos $ srn $13 — 
—0.1 and Ax 2 — 0.2 at $13 = (see the discussion in Ref. [31]). The contrast is due to 
the different marginalization procedures in the frequentist and Bayesian theories: in the 
frequentist theory only the minimum of x 2 with respect to the marginalized parameters is 
considered, whereas in the Bayesian theory marginalization is implemented by integrating 
over the distribution of the marginalized parameters, as we have done, for example, in 
Eq. (13. 3p . In the case of the marginalization over cos 5 = ±1, the Bayesian procedure of 
summing the prior probability distribution over cos 5 = ±1 in Eq. (14.11) for each value of 
sin?!^ is different from the frequentist consideration of ^(sin^) f° r cos 5 = —1 only, 
which is due to x 2 ( srn "$13 1 cos 5 = — 1) < x 2 ( sm ^i3| cos 5 = 1). In general, the Bayesian 
marginalization procedure has the merit to take into account all the distribution of the 
marginalized parameters, which gives more information than the single point of minimum 
of re- 
using the informative prior on $ 13 in Eq. (14.11) . from the analysis of solar and Kam- 
LAND data We found the best-fit point, corresponding to the maximum of the posterior 
probability distribution, 

Ar/4 = 7.58 x 10~ 5 eV 2 , sin 2 tf 12 = 0.31 , sin 2 tf 13 = 0.012 . (4.3) 

The shadowed areas in Figs. [TOl and [TT1 show the posterior credible regions with 90%, 95%, 
and 99.73% probability in the sin 2 ^^-Am^ and sin 2 $ 12 -sin 2 planes, respectively. 
The boundaries of the corresponding regions obtained with an uninformative prior, shown 
in Figs. [5] and El are depicted with solid lines. 

Since the prior information constrains only $13, the best-fit values and allowed ranges 
of Am 2 and $12 are similar to those obtained in Section [3] with an uninformative prior. 
A small change is due to the correlation with $ 13 . 
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Figure 10: Shadowed areas: 90%, 95%, and 99.73% Bayesian credible regions in 
the sin 2 , &i 3 -Am^ 1 plane obtained with the informative prior probability distribution in 
Eq. (14. ip . The regions enclosed by solid lines correspond to those in Fig. 0, obtained with 
an uninformative prior. The straight vertical dotted lines enclose, respectively, the 90%, 
95%, and 99.73% prior credible regions of sin 2 $13. The left and right plots correspond, 
respectively, to a normal and an inverted scheme (see Fig. [1]). The cross and asterisk indi- 
cate, respectively, the best-fit points of the analyses with uninformative and informative 
priors. 



On the other hand, one can see that the assumption of the informative prior in 
Eq. (14. ID leads to a significant reduction of the allowed range of $13 with respect to that 
obtained with an uninformative prior, as should have been expected. 

A curious feature of Figs. Hn]and[TT]is that the 90% and 95% allowed ranges of sin 2 $ 13 
seem to be larger than those allowed by the prior distribution (vertical straight lines in 
Fig. [inland horizontal straight lines in Fig. fl~T|) . Such a conclusion would be erroneous, 
because the prior distribution in Eq. ( 14. ip concerns only one parameter, whereas the 
credible regions in Figs. [10] and [11] constrain two parameters taking into account their 
correlation. 

The posterior probability distribution of sin 2 $ 13 obtained from the marginalization 
in Eq. ( 13. 3p implies an allowed range of sin 2 $13 which is smaller than that given by the 
prior distribution, as one can see from Fig. [12] We obtained the 90%, 95% and 99.73% 
upper bounds 

. 2 f 0.027, 0.030, 0.045 (normal scheme) , . . 

sin it 13 < I Q 03Q ^ Q Q33 ^ Q 04g (i nver ted scheme) , ^ ' 

which are smaller than the corresponding ones in Eq. ( 14. 2ft . These bounds are also about 
60% smaller than those obtained in Eq. (13 ,4p with an uninformative prior. Figure H3] 
shows the comparison of the posterior probability with the one in Fig. [HI which has been 
obtained with an uninformative flat prior. 

It is interesting to note that the bounds on sin 2 $13 in Eq. ( I4.4p are similar to those 
obtained with a global x 2 analysis of neutrino oscillation data in Ref. [9] (see, however, 
the caveat on the comparison of frequentist and Bayesian results discussed at the end of 
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Figure 11: Shadowed areas: 90%, 95%, and 99.73% Bayesian credible regions in the 
sin 2 ^ 12 -sin 2 $ 13 plane obtained with the informative prior probability distribution in 
Eq. (14. ip . The regions enclosed by solid lines correspond to those in Fig. El obtained 
with an uninformative prior. The straight horizontal dotted lines enclose, respectively, 
the 90%, 95%, and 99.73% prior credible regions of sin 2 $13. The left and right plots 
correspond, respectively, to a normal and an inverted scheme (see Fig. CD). The cross and 
asterisk indicate, respectively, the best-fit points of the analyses with uninformative and 
informative priors. 

Section [3]). Our results also agree with the weakening of the hint of $13 > discussed 
in Ref. [9] coming from the addition of atmospheric and long-baseline accelerator and 
reactor neutrino data to the analysis of solar and KamLAND data: using the method 
described in Section [31 the significance of the hint of $13 > is reduced from about 1.2a to 
about 0.8cr (with 0.72 and 0.75 respective probabilities of the smallest posterior credible 
region which includes $13 = in the normal and inverted schemes). The discrepancy 
with the 1.6cr reported in Ref. [27] is probably due to the marginalization over cos 5 = ±1 
in Eq. (14.11) . In fact, for cos 5 = —1 Fig. 24 of Ref. [31] implies a prior in favor of -#13 > 0, 
which leads to a global hint of $13 > at the 1.5cx level (with 0.93 probability of the 
smallest posterior credible region which includes $13 = 0), in agreement with the value 
in Ref. [27] (1.6<r). Let us however emphasize that, since the marginalization over the 
unknown value of 5 is the correct procedure in the Bayesian approach, our result for the 
statistical significance of the global hint of $13 > is 0.8a. 

Let us finally remark that the results presented in this Section depend on the choice of 
the prior probability distribution for $i 3 obtained from the fit of the data of atmospheric 
and long-baseline accelerator and reactor neutrino experiments. In Eq. (14. ip . instead of 
the x 2 °f Ref- [31] we could have used, for example, the x 2 °f one °f Refs. [9,32-34]. 
However, since in these papers the same data have been fitted with similar assumptions 
and methods, using one of these x 2 's would not change dramatically the numerical results 
presented above. For example, we considered the x 2 m Fig- 3 of Ref. [9], which corresponds 
to the prior upper bounds 

sin 2 tfu < 0.030 (90%) , 0.039 (95%) , 0.063 (99.73%) . (4.5) 
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Figure 12: Marginal posterior probability distribution of sin 2 $13 obtained with the 
informative prior probability distribution in Eq. (14.11) . The dotted curve shows the prior 
distribution in Eq. (14. ip . The long (short) straight vertical lines show the 90%, 95%, and 
99.73% posterior (prior) probability levels. The short straight dotted vertical lines have 
been slightly shifted to the right to avoid superposition (compare with Eqs. (14.21) and 

D)). ' 

We obtained the best-fit values 

Am 2 21 = 7.58 x l(T 5 eV 2 , sin 2 t? 12 = 0.31 , sin 2 t? 13 = 0.008 , (4.6) 

and the posterior upper limits 

sin 2 tf 13 < 0.030 (90%) , 0.033 (95%) , 0.051 (99.73%) . (4.7) 

One can see that these values are close to the corresponding ones in Eqs. f l4.2l) - fl4.4|) . 
For the hint of $13 > we have a 0.9a statistical significance (with 0.78 probability of 
the smallest posterior credible region which includes $13 = 0), in perfect agreement with 
Ref. [9]. 

5 Conclusions 

In this paper we presented the results of a Bayesian analysis of the solar and KamLAND 
neutrino data with the aim of determining the value of the unknown mixing angle $13 in 
the framework of three-neutrino mixing. 

We found that with an uninformative flat prior distribution in the relevant mixing 
parameters Am 2 2 , sin 2 -^, sin 2 ^, the Bayesian credible regions in the sin 2 $ 13 — Am 2 2 
and sin 2 $12 — sin 2 $13 planes are only slightly smaller than the allowed regions obtained 
with a traditional least-squares analysis, implying a rather stringent upper bound for 
sin 2 $ 13 . Our analysis confirms the 1.2cr hint of $13 > found in Ref. [27]. 

We also performed an analysis with an informative prior which represents informa- 
tion on ^ 13 obtained in atmospheric and long-baseline accelerator and reactor neutrino 
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Figure 13: Marginal posterior probability distribution of sin 2 $13 obtained with the infor- 
mative prior probability distribution in Eq. (14.11) . The dotted curve shows the marginal 
posterior distribution in Fig. [HI obtained with an uninformative flat prior. The long 
straight vertical lines show the 90%, 95%, and 99.73% posterior probability levels in 
Eq. (14.41) . The short straight vertical lines show the corresponding probability levels in 
Eq. (13.41) . obtained with an uninformative flat prior. The short straight dotted vertical 
line in the figure on the right has been slightly shifted to the right to avoid superposition 
(compare with Eqs. (13 .4p and (14.41) ). 

experiments, independently from solar and KamLAND neutrino data. We found that 
such a prior implies a significant decrease of the upper bound on $13 with respect to that 
obtained with an uninformative prior and the hint of $13 > is reduced to a 0.8<j level. 
Our results are similar to those obtained with a global \ 2 analysis of neutrino oscillation 
data in Ref. [9] (see, however, the caveat on the comparison of frequentist and Bayesian 
results discussed at the end of Section [3]). 

Let us finally emphasize that Bayesian inference (see Refs. [35-39]) is founded on a 
consistent theory and can always be implemented in a correct way (given enough com- 
putational power). On the other hand, the frequentist method is based on an unphysical 
definition of probability and in most cases of interest cannot be implemented in a correct 
way. In particular, a dramatic flaw of the frequentist method is that the frequentist def- 
inition of probability does not allow the treatment of theoretical and systematic errors 
as random variables. Hence, the aim of the frequentist statistics approach of extract- 
ing objective statistical information from data cannot be realized in practice. Since the 
Bayesian theory does not suffer from such shortcomings, we think that it is preferable for 
attaining reliable results from the analysis of experimental data. 

A Regeneration of Solar z/ e 's in the Earth 

Solar neutrinos arriving at a detector during night-time pass through the Earth, where 
the matter effect (also called "MSW effect" [40,41]) can cause a change in the flavor 
composition, which is called "regeneration of solar i/ e 's in the Earth". In this Appendix, 
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we derive the connection between the averaged probability of survival of solar electron 
neutrinos passing through the Earth, P 1 ^L fV , which is measured during night-time, 
the averaged probability of u e survival from the core of the Sun to the surface of the Earth, 
Pj^v i which is measured during day-time, and the probability of v 2 — > v e transitions in 
the Earth, P^ 8 ^. We also discuss the connection between P^£ in the case of three- 
neutrino mixing and the probability of v 2 — > v e transitions in the Earth in the case of 
two-neutrino mixing. 

The mixing of neutrino states is given by 

3 

K) = ^2 U akWk) (a = e,n,r), (A.l) 

k=l 

where U is the 3x3 unitary mixing matrix of the neutrino fields (see Refs. [1-3]). We 
adopt the standard parameterization [4, 5] 

C12C13 S12C13 si 3 e~" A 

U = I -S12C23 - ci 2 s 2 3S 13 e t5 c 12 c 23 - si 2 s 2 3Si 3 e t5 S23C13 I , (A.2) 

■512-523 ~~ Ci 2 C 23 Si 3 e tS — C12S23 — Si 2 C 2 3Si3e l5 C23C13 

where c a b = cos$ a b and s a b = sin^. The three mixing angles $12, $13, $23 take values 
in the ranges < d a b < tt/2. The CP- violating phase 5 is confined in the interval 
< 5 < 2tt. We neglected possible Majorana phases, which are irrelevant for neutrino 
oscillations [42-44]. 

A solar neutrino, created in the core of the Sun as a u e , arrives in a detector as a 
superposition of Vi, u 2 , and z/3. However, since the neutrino squared-mass differences are 
relatively large (see Eqs. fll.ll) - fll.5l) ). the average of the oscillation probability over the 
energy resolution of the detector washes out the interference terms between the massive 
neutrinos [45]. This is due to the fact that the vacuum oscillation lengths are much 
shorter than the Sun-Earth distance: 

AttE ( E \ AttE ( E \ 

L°l c = — -^~30km — — , L™ c ~ L°l c = — ^~lkm — — , (A.3) 

where E is the neutrino energy, which in solar neutrino experiments varies in the interval 

0.2 MeV < E < 15 MeV . (A.4) 

Then, the measurable averaged survival probability of solar electron neutrinos after cross- 
ing the Earth is given by 



-pSun+Earth ^ Sun p E arth /a c\ 

k=l 



where P^+ v . is the probability of v e — > transitions from the solar core to the surface 
of the Earth and P^ff^ e is the probability of — > v e transitions in the passage through 
the Earth. 
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In matter, electron neutrinos feel a charged-current potential Vcc — V^GpN e , where 
Gp is the Fermi constant and N e is the electron number density. The quantity which 
gives the matter effect in the evolution equation of neutrino flavors is 



Acr, = 2EV, 



cc 



1.53 x l(T 7 eV 



N\ cm - 



(—) 

VMeVy 



(A.6) 



where Na is the Avogadro number. The electron number density in the solar core is about 
lOOA^Acm -3 . In the Earth, the electron number density varies from about 2.2 Aachi -3 
in the mantle to about 5.5 A^a cm -3 in the core. Thus, for solar neutrinos we have Aqc ^ 
2.3 x 10 _4 eV 2 , which is much smaller than the atmospheric squared-mass difference 
Am ATM (see Eq. ( jl.2p ). This means that the matter effect cannot induce transitions 
between z/ 3 and the two neutrinos V\ and since the two groups are separated by the 
large atmospheric squared-mass difference Am ATM (see Eq. (11.51) and Fig. CE]). In other 
words, the massive neutrino component z/ 3 propagates without disturbance from the core 
of the Sun to the detector and the corresponding transition probabilities in Eq. (1A.5j) are 
simply given by 

(A.7) 



P Sun^ =p ™ =1(^)^=1^1 



jEarth 



Furthermore, taking into account the conservation of probability, we have 



^Sun 

v e — YV\ 
pEarth 

Vl~*V e 



P, 



Sun 

Ve~*V Z 

Earth 



Pi 

R 



Sun 
Earth 

V 2 ~*Ve 



\u*\ 
\u*\ 



pSun 

pEarth 
*V2-*Ve ' 



(A.8) 
(A.9) 



•Sun 



Let us now express the averaged survival probability P v ^ v of electron neutrinos from 



the solar core to the surface of the Earth in terms of the transition probabilities P^ 



un 
v e — >Vh ' 



-^Sun 

* v e — »f e 



l<fe|<Sk>| 2 



y^ y (VeWk)(Vk\SWe) 
k=l 



(A.10) 



k=l 



where S is the evolution operator. We neglected the interference terms for the reason 
discussed above, before Eq. (IA.5I) . Using Eqs. (IA.7I) . flA.8j) . and flA.101) . we can express 



P^ u H, Jn in terms of P ™. : 

V e — >V2 V e —¥V e 



;>Sun 

V e — *V 2 



\U el \ 2 (l-\U e3 \^ 



\Ue^~P 



■Sun 



\Uel\ 2 -\U e 2\ 2 



(ah; 



Finally, using Eqs. (tOl . ([Oil . itOl) . and (lATTT) . we obtain, from Eq. (fA~5l) 



-=rSun+Earth -=rSun 

P v e — >V e Ve— >V e ~^ 



a- 


u e3 


2 f - 2 (P^Ve ~ 


u e3 


') 




I" pEarth 

L V2— >V e 


Ue2\ 2 } 




u el \ 2 - 


u e2 


2 



(A.12) 



Since in practice \U e i\ 2 > \U e2 \ 2 , because sin 2 $ 12 < 1 (see Refs. [10,31]), and \U e3 \ 2 
is small, there is a regeneration of electron neutrinos in the Earth if PR 8 ^ > \U e 2\ 2 . 
Note that in the absence of matter effects, we have P^f^ = K^l^e)] 2 = \U e 2\ 2 and 

-=rSun+Earth -=rSun 

V e ^l>e Ve— *V e ' 
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Let us now discuss the calculation of P,^ 3 ^ . The evolution of neutrino flavors in 
matter is governed by the Schrodinger equation (see Refs. [1-3]) 

d 

i — — \l/p = Hp \l/p , (A.13) 

ax 

with the effective Hamiltonian 

Hp = — ([/AM 2 W + A) , (A. 14) 

and 

/0 \ [A cc 0\ 

AM 2 = Am 2 ! , A= 0. (A.15) 
\0 Am^/ \ 0/ 

Here, ^ Q = (^ a |z/) is the amplitude of the flavor a in the state \v) which describes 
a propagating neutrino. The column matrix ^/p of flavor amplitudes is related to the 
column matrix \I>m = (^i, i>2, ipz) T of mass amplitudes (ipk = (vk\v)) by 

y F = U^ M . (A.16) 

In the calculation of P^ 3 ^, the initial mass and flavor amplitudes are 

<M0)=5 fc2 , M0) = U a2 - (A.17) 

The probability of z/ 2 — * transitions at a distance £ from neutrino production is given 
by 

Pv n (x) = \M%)\ 2 ■ (A.18) 

Taking into account the fact that the mixing matrix in the parameterization in 
Eq. ( 1A.2I) can be written as 

U = R 23 W 13 R 12 , (A.19) 

with 

/ ci 2 s 12 0\ /l \ / ci3 s 13 e- i>N 

P 12 = -s 12 c 12 , P 23 = c 23 s 23 , W 13 = 1 
\ 1/ \0 -s 23 c 23 / \s 13 e l5 cia 

(A.20) 

it is convenient to work with the new column matrix of amplitudes \I/ = (^1,^2, ^s) T 
defined by 

$ = p^ 13t P 23t ^ F = P 12 ^ M , (A.21) 
which follows the evolution equation 

iA$ = H$. (A.22) 
ax 

Since R 23 commutes with the matter potential matrix A, the new effective Hamiltonian 
EI is given by 

§ = — (P 12 AM 2 P 12t + W l3] A W 13 ^j . (A.23) 
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Explicitly, we have 



1 



's?2 Am 2i + c? 3 ^cc Ci 2 Si 2 Am^ -c 13 si 3 e l& A cc 



H = — c 12 s 12 Am 2 21 c 2 12 Am 2 21 . (A.24) 



2E 



-ci 3 si 3 e iS A C c Amt + s 2 3 Acc 



From Eq. (1A.21[) . we have ip 3 = ip 3 . Therefore, ip 3 is the amplitude of v 3 . Since 
Am| x ^> Acc-, in practice the third eigenvalue of EI is equal to Am 2 1 /2E and the mat- 
ter effect cannot induce transitions between z/ 3 and the other two massive neutrinos, as 
discussed above. Furthermore, since ^(0) — ^3(0) = (from Eq. (1A.17I) ). in practice 
the contribution of u 3 is negligible and P^f^ can be calculated by solving the effective 
two-neutrino evolution equation 

i-^-^ = m, (A.25) 

da; 

with \1> = {ipi,ip 2 ) T = {ipi,ip2) T and 

1 ( s\ 2 Am 2 2l + c 2 13 A cc Ci 2 s 12 Am 2 2l 



m = — ^ 2 """ 21 . ^ 2 a 2 • ( A - 26 ) 

2E V c 12 Si 2 Am 2 2l c\ 2 Am 2 2l J v ; 

This effective Hamiltonian coincides with the effective Hamiltonian in the case of two- 
neutrino mixing (see Refs. [1-3]), with the matter contribution Acc multiplied by the 
three-neutrino mixing factor cf 3 . The initial column matrix of amplitudes is explicitly 
given, from Eqs. (1A.17I) and (1A.21I) . by 

*<°> = {z :::) (!) = (:::) • ^ 

The probability of v 2 — > v e transitions at a distance x from neutrino production is given 

P^(x) = \[R 23 W 13 $(x)} e \ 2 = c\,mx)\ 2 . (A.28) 
Therefore, in practice, the probability of v 2 — > u e transitions in the Earth is given by 

P££. = (1 - l^sl 2 ) , (A.29) 

where P^^ 2v is the probability of v 2 — > u e transitions calculated in the case of two- 
neutrino mixing with an effective matter contribution multiplied by c 2 3 = 1 — |t/ e3 | 2 . 
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